Magnetic susceptibility of Cerium: an LDA+DMFT study 
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The magnetic properties of Ce in the a and 7 pliase are calculated within the LDA+DMFT 
approach. The magnetic susceptibility in these two phases shows a similar behavior over a wide 
temperature range: a Curie- Weiss law at high temperatures, indicating the presence of local mo- 
ments, followed by a maximum in a crossover regime, and a saturation characteristic of a state with 
screened local moments at low temperature. The difference in experimentally observable magnetic 
properties is caused by the shift of the susceptibility to higher temperatures in the Q-phase compared 
to the 7-phase. 

PACS numbers: 71.27.-|~a, 71.20.Eh 
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I. INTRODUCTION 

The isostructural a — 7 transition in Ce is one of the 
classical problems in modern solid states physics. In the 
low temperature a-phase (up to T~ 100 K at normal con- 
ditions, or up to T~ 300 K for a pressure P of 1 GPa) Ce 
behaves like a Pauli paramagnet, while in the high tem- 
perature 7-phase the susceptibility approximately follows 
a Curie- Weiss lawii The transition is accompanied by a 
drastic volume collapse (9-15%)^ and dramatic changes 
of the electronic spectrai^ 

A number of theoretical models were proposed to de- 
scribe the q: — 7 transition. One of the first was a promo- 
tional model, where localized 4/ electrons were suggested 
to transfer to the spd— (valence) band state, losing their 
local momentsj^i^ This was in contradiction to later ex- 
perimental results that showed that the number of 4/ 
electrons is almost unchanged during the transitioui^ As 
a result, a Mott-like picture was proposed, where the va- 
lence of the Ce ions does not change, but the transition, 
which affects the degree of 4/ electron localization, oc- 
curs as a result of the change of the ratio of on-site / — / 
Coulomb interaction (U) to kinetic energy<^ 

Further neutron experiments^ confirmed that the Ce- 
4/ electrons remain localized and indicated that the 
Kondo volume collapse model^ii^ could be more plausi- 
ble. According to this model the 4/ electrons remain lo- 
calized in the low temperature a— phase, but form Kondo 
singlets with conduction electrons in the spd— hand. As 
a result of this strong coupling the local spin moments of 
the 4/ electrons get screened, which leads to the volume 
collapse in the Kondo regime. 

Thus a full account of the hybridization between lo- 
calized 4/ and band spd electrons is need for the cor- 
rect description of Ce. Moreover, it was recently shown 
that / — / hopping is also important and should be taken 
into consideration The full information about the non- 
interacting band structure of Ce can be obtained in the 
frameworks of the density function theory (DFT), e.g. 
in the Local density approximation (LDA). These den- 



sity functional calculations can be extended to include 
local correlations within the LDA-f-DMFT scheme (com- 
bination of local density approximation and dynamical 
mean- field theory) 

"Ab-initio" calculations in the LDA-I-DMFT approach 
were successfully applied to the modeling of electronic 
and structural properties of Ce. They have illustrated 
the key role of the entropic contribution to the free cn- 
g,j.g yi2^i3 g^j^j |-]^g iniportance of the formation of a quasi- 
particle peaki^i^ for the description of the a — 7 tran- 
sition in Ce. The applicability of the Kondo model was 
mainly discussed via an analysis of the spectral proper- 
ties of the two phases: the temperature dependence of 
the quasi-particle resonanc o^^'^^ or features of the Ce- 
spd bandsii. Meanwhile, the description of a key phys- 
ical observable, namely the temperature dependence of 
magnetic susceptibility (for which different models were 
originally proposed) has not yet been attempted. 

In the present report we use the LDA-t-DMFT method 
to calculate the magnetic properties of Ce in the a and 7 
phases. The results of the study show that the temper- 
ature dependence of the magnetic susceptibility is very 
similar in Ce— a and Ce— 7 and that both phases should 
thus be described by the same model. The difference in 
the observed magnetic properties is attributed to the de- 
crease of the hybridization between localized / and con- 
ductive spd electrons, which leads to a shift of the mag- 
netic susceptibility curve to lower temperatures in Ce— 7. 



II. CALCULATION DETAILS 

We performed LDA calculations using the Linearized 
muffin-tin orbitals (LMTO) method An almost orthog- 
onalized version of the LMTO in the F- representation 
with Ce - 6s, 6p, 5d, and 4f states included to the basis 
set was used. The Hamiltonian was generated on a mesh 
of 1728 fc— points in the full Brillouin zone (BZ). This 
LDA hamiltonian was then transformed to a basis set 
with a diagonal form at the F-point. In this basis set the 
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FIG. 1: Uniform magnetic susceptibility xiT) for Ce—a. In- 
set (a) - enlarged view of x{T) in the low temperature region. 
Inset (b) - inverse magnetic uniform susceptibility 



three lowest energy states at the F point correspond to 
the tiu, the next three to the <2ti and the highest energy 
states to the 0211 irreducible representation. 

The on-site Coulomb repulsion parameter (U) was es- 
timated to be 6.0 eV using a constrained supercell calcu- 
lation. This is in agreement with previous resultsJ^ The 
intra-atomic Hund's rule coupling was set to Jh=0 eV. 

For the solution of the DMFT equations we employed 
a diagrammatic ('continuous-time CT-HYB') quantum 
Monte Carlo algorithm which samples the partition func- 
tion in powers of the impurity-bath hybridization] ^"-"i^ ^ 
The Coulomb term was treated in the density-density 
form. The double counting correction was set to E^c = 
[/(riDMFT ~ 5))^ with TiDMFT the total number of 4/ 
electrons self-consistently obtained within DMFT. The 
LDA-I-DMFT calculations was not fully self-consistent in 
the sense that the LDA charge density was not recalcu- 
lated after the DMFT run. This can be done since the 
total number of the 4/-electrons doesn't change signifi- 
cantly. The Ce— 4/ spectral functions were calculated us- 
ing the maximum entropy methodi^ We first computed 
the self-energy S on the real frequency axis and used it to 
obtain the orbitally resolved and total spectral functions. 

We calculated the uniform magnetic susceptibility x as 
the ratio of the field-induced magnetization dm(T) and 
the energy change SE associated with the applied field 
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Here n-^ and nj, are the total occupation numbers for spin 
up and down. The susceptibility x '^as calculated for 
SE = 0.01 eV, which is within the interval where ■m{T) 
is a linear function of SE. 



FIG. 2: (Color online) Imaginary part of the self-energy for 
all 14 4/-orbitals at T=129 K (/3 = 90 eV"^) in Ce-a. Due 
to the crystal-field splitting they form three different sets of 
curves denoted as tiu, t2u and a2u- The inset shows the low- 
energy behavior of ImE(ia'„). For one of the curves exact 
positions of Matsubara frequencies and linear extrapolation 
(by dashed line) to zero are shown. 



III. RESULTS 



The uniform magnetic susceptibility x(2^) ^oi the Ce-a 
phase, obtained using the LDA-I-DMFT calculations, is 
shown in Fig. [1] In the low temperature region the mag- 
netic susceptibility is temperature independent up to ~ 
300 K. This is in qualitative agreement with experimental 
findingsi^ The absolute value of the magnetic suscepti- 
bility is underestimated in LDA-I-DMFT by comparison 
with X measured for P=l GPa, where the presence of the 
structurally different Cc-/3 phase is minimal: The experi- 
mental value of X is 0.6 — 0.7 x 10^'^ cmu/mol, while the 
theoretical estimate is 0.22 x 10^"^ emu/mol. The under- 
estimation of the magnetic susceptibility could be due to 
the absence of long range correlations in our single-site 
DMFT calculations. Also, the spin-orbit coupling and 
associated orbital moment are not captured within our 
simple LDA description. 

The plateau in the susceptibility of Ce at low- 
temperatures corresponds to a coherent regime, where 
all 4/ electrons with local moments are screened. In- 
deed, one may notice from Fig. [2] that in this tempera- 
ture region the imaginary part of the self-energy for all 
4/ orbitals is linear and approaches zero for low Matsub- 
ara frequencies. This is a signature that the system at 
these temperatures is in a coherent Fermi-liquid regime. 

The importance of the screening effects by spd elec- 
trons is most clearly seen when one compares Fig. [U 
inset (a) and Fig. [3] The local magnetic susceptibility 
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FIG. 3: Local magnetic susceptibility xo{T) and inverse mag- 
netic local susceptibility Xo^C^) Ce-a. 

presented in Fig. [3] was calculated as 

Xo= f\s,{T)SM)dr, (2) 
Jo 

where /3 is the inverse temperature, and {Sz{t)Sz{0)) is 
the imaginary time dependent spin-spin correlation func- 
tion. In spite of the strong non-linearity of the inverse 
uniform susceptibility, the local susceptibility follows a 
Xo^{T) = T/C law up to T« 1000 K. This implies that 
the 4/-electrons retain the local nature of the magnetic 
moment even at very low temperatures and there is no 
dclocalization process. The violation of the Curie- Weiss 
law of the uniform susceptibility (Fig. [l] inset (b)) and 
coherence of the system (which is seen in Fig.[2|) is caused 
by the screening of local spins by conduction electrons. 

With increasing temperature the screening effects 
weaken. This leads to an increase of the entropy accord- 
ing to Ref. [l^. Note that the growth of the susceptibility 
starts at 350 K, i.e. in the region where the a — 7 
transition occurs. 

At temperatures above 350 K, the experimentally sta- 
ble phase of Ce is the 7, not the a phase. However, we 
may still simulate the magnetic properties of Ce— a even 
in a temperature region, where it does not exist by us- 
ing the LDA band structure of the corresponding phase. 
Making such calculations one concludes that the tem- 
perature range 350 K<T<1000 K is a crossover region, 
where the system can neither be described as a coherent 
electron liquid nor as a lattice of localized /-states weakly 
hybridized to conduction electrons. The maximum in the 
uniform magnetic susceptibility at 820 K appears as a 
result of a competition between Kondo and local spin 
regimes. Interestingly, in the Coqblin-Schricffer model 
this maximum appears if the degeneracy of the impu- 
rity level is more than 3^^^ This is in agreement with 
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FIG. 4: Uniform magnetic susceptibility xi^) and inverse 
magnetic uniform susceptibility x~^{T) for Ce-7, inset (a). 
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FIG. 5: Local magnetic susceptibility Xo{T) and inverse mag- 
netic local susceptibility Xo^C^) ^'^^ Ce-7. 



the present results, where the degeneracy of the lowest 
energy <i„ states equals 6 (for both spins). 

A behavior consistent with the Curie- Weiss law, 
l/x(T) ~ T, for the uniform susceptibility is only seen 
above 1000 K (Fig.[Tl inset (a)). 

The uniform magnetic susceptibility in the 7 phase is 
qualitatively similar to the one just described for Cc— a. 
A Curie- Weiss like behavior at high temperature, then a 
crossover region characterized by a maximum in the sus- 
ceptibility, followed by a drastic drop and (at the lowest 
accessible temperatures) a saturation. The difference is 
mainly in the numbers. Already at T=300 K, the uniform 
magnetic susceptibility of Ce in the 7-phase is Curie- 
Weiss like. The broad maximum marking the crossover 
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FIG. 6: (Color online) Imaginary part of the self-energy for all 
14 4/— orbitals at T=64 K in Ce— 7. Due to the crystal-field 
splitting they form three different sets of curves denoted as 
till, t2u and a2u- The inset shows the low-energy behavior of 
luiE{iijj„). For one of the curves exact positions of Matsubara 
frequencies and linear extrapolation (by dashed line) to zero 
are shown. 



region is at ^ 200 K, much lower than in Ce— a. How- 
ever, since the 7— phase experimentally exists only above 
300 K the crossover regime and the constant x(T) regime 
found below ~90 K are not accessible in measurements. 

The absolute value of the calculated x in the region 
around 500 K, where it is experimentally measurable, is 
underestimated (like in the a— phase). We find ^ 0.7 x 
lO^'^ cmu/mol against the experimental ^ 1.4 x 10""^ 
emu/mol2^. As one may see in Fig. [5] the local magnetic 
susceptibility, xoj for Ce— 7 is similar to that for Ce— a 
and follows a C/T law up to T < 1000 K. 

At higher temperatures Xo^ shows a deviation from 
the linear behavior with an upturn of the curve. A simi- 
lar upturn can be observed for the inverse susceptibility 
of the one-band Hubbard model in the correlated metal 
regime (U=2.5, W=2)^ This effect is due to the fact 
that at higher temperatures the thermal fluctuations be- 
come so large that the spins at different r points be- 
come uncorrelated and "do not feel" each other, i.e. that 
(S'^(t-)^^(O)) - C with T ^ 00. The resuh is a 
deviation of the local susceptibility Xo^ = T/C{T) from 
the Curie law. 



IV. DISCUSSION 

The shape of the susceptibility curves for the a- and 
7-phases is very similar. Moreover, it strongly resem- 
bles x{T) obtained in models with a projected Hilbcrt 
space, such as the Coqblin-Schricffcr impurity ;^"^'^^ its 
lattice version^ and also x(T) in the Periodic Anderson 
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FIG. 7: (Color online) Top panel: Spectral function for the 
Ce— a phase calculated for P — 90 eV~^. Bottom panel: Spec- 
tral function for the Ce— 7 phase calculated for /3 = 180 eV~^ 



modeli^ The uniform static magnetic susceptibility in 
all these models is characterized by a high-temperature 
Curie- Weiss tail, a maximum in the intermediate regime 
and then a drastic drop with a constant susceptibility at 
the lowest temperatures. 

The high-temperature behavior of xi'^) is explained by 
the presence of local moments and described in a similar 
manner in all of the models, while the drastic decrease of 
the magnetic susceptibility at intermediate temperatures 
is caused by the screening of this moment. The lattice 
effects such as non-zero hopping between localized and 
conductive states, centered at different sites are already 
taken into account in our LDA-I-DMFT calculation. The 
results clearly show the formation of a coherent state with 
Iml](ia;„) — >■ at low temperatures (see Fig. l2l6p . Mean- 
while, an inspection of the spectral functions, plotted in 
Fig. [71 shows that there is no gap or pseudo-gap in the 
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vicinity of the Fermi level in the /-electron spectral func- 
tion. Note that we present here results at considerably 
lower temperatures than in previous DMFT studies. 

Comparing Fig. [T] and H] one sees that the behavior of 
the uniform magnetic susceptibility in a— and 7— Ce is 
qualitatively the same. The susceptibility in the 7 phase 
seems to be shifted to the low-temperature region and 
rcnormalizcd (as compared with Cc— a). Thus, one may 
argue that those phases are physically similar in a wide 
temperature range. 

The change of the lattice volume under the a — 7 tran- 
sition results in a modification of the hybridization func- 
tion in the vicinity of the Fermi level and decrease of the 
/ — / hopping. The weakening of the hybridization in 
Ce— 7 (see Ref. [iBfor instance) leads to a decrease of the 
exchange parameter Jk in the Kondo model according 
to the Schrieffer- Wolff transformation.— Numerical cal- 
culations show that the decrease of the Kondo exchange 
results in a shift of the maximum of the magnetic suscep- 
tibility to lower temperaturesi^ This is exactly what is 
seen in Ce under the a — 7 transition. A similar shift of 
the maximum of x(^) due to a change of the hybridiza- 
tion was also observed in the study of the magnetic prop- 
erties of the Periodic Andersen model (PAM)<2^ 

It is also instructive to compare the present results 
with the situation in Pu, which is on the border of 
the transition between the actinide elements with local- 
ized and delocalized electronsi^ Moreover the degree of 
localization changes in the different phases of Pu<^i^ 
The DMFT calculations of Pu show that the local mag- 
netic susceptibility dramatically changes from Pauli-like 
to Curie- Weiss as the volume increasesi^ This is com- 



pletely different from the situation observed in Ce, where 
the local susceptibility follows a 1/T law both for the a 
and 7 phases. This is further evidence for the absence of 
any localization/delocalization transitions in the 4/-shell 
of Ce, i.e. against the Mott transition scenario. 

To summarize, in the present paper the magnetic sus- 
ceptibility for the a and 7 phases of Ce was investigated. 
It exhibits a qualitatively similar behavior in both phases. 
There is no Mott transition for any of the Ce phases. On 
the local level they are both characterized by the presence 
of magnetic moments, screened by band spd electrons 
at low temperatures. With increase of the temperature 
this coherent state with constant susceptibility is grad- 
ually destroyed by thermal excitations, which results in 
the formation of a Curie- Weiss paramagnetic state in the 
high temperature region. The difference in experimen- 
tally observable magnetic properties for the two phases 
of Ce is related to a shift of the susceptibility in Ce— 7 
to lower temperature region (as compared to a phase). 
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